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Abstract. An analysis of the high-energy emission from IGR J16393-4643 (=AX J 1639.0-4642) is presented using data from 
INTEGRAL and XMM-Newton. The source is persistent in the 20-40 keV band at an average flux of 5.1 x 10~"ergs cirT 2 s , 
with variations in intensity by at least an order of magnitude. A pulse period of 912.0±0.1 s was discovered in the ISGRI and 
EPIC light curves. The source spectrum is a strongly-absorbed (A'h = (2.5 ± 0.2) x 10 23 cnr 2 ) power law that features a high- 
energy cutoff above lOkeV. Two iron emission lines at 6.4 and 7.1 keV, an iron absorption edge >7.1 keV, and a soft excess 
emission of 7 x 10~ 15 ergs cnr 2 s -1 between 0.5-2 keV, are detected in the EPIC spectrum. The shape of the spectrum does not 
change with the pulse. Its persistence, pulsation, and spectrum place IGR J16393-4643 among the class of heavily-absorbed 
HMXBs. The improved position from EPIC is RA. (J2000)= 16 h 39 ra 05.4 s and Dec.= -46°42'12" (4" uncertainty) which is 
compatible with that of 2MASS J16390535-4642137. 
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1. Introduction 

! The INTEGRAL core program (CP: IWinkler et alll2003h rou- 
tinely devotes observation time to Galactic Plane Scans (GPS) 
and Galactic Centre Deep Exposures (GCDE). These numer- 

" ous snapshots of the Milky Way can be assembled into mo- 
saic images of long exposure time (~1 M s). This gamma-ray 
view of the galaxy, as collect ed by ISG RI ( Uber tini et all2 003: 
lLebrunetalJl2003l) . enabled iBird et al J J2004 to detect 123 
sources at a significance above 6cr. Around 20 of these sources 
are of unknown origin. A good portion of these new sources 
probably belong to the class of heavily-absorbed High Mass X- 
ray Binaries (HMXBs) that are concentrated along the galactic 
plane and in the spiral arms. 

High-Mass X-ray Binaries are composed of a compact ob- 
ject such as a neutron star or a black hole that orbits a mas- 
sive stellar companion. Depending on the type of comp anion, 
know n HMXBs can be divided int o two groups Jvan Par adiis 
Il983l) . Most HMXBs classified bv lLiuet alJ d200(j) contain a 
Be star. These systems are usually transient sources with hard 
spectra. The compact object has a wide orbit which mostly 
keeps it away from the Be star and its disk. Outbursts in these 
systems are due primarily to the compact object approaching 

Send offprint requests to: arash.bodaghee@obs.unige.ch 



the star and accreting matter from the slow, dense stellar wind. 
The second group of HMXBs features an O or B supergiant 
star. The orbit of the compact object places it well within the 
stellar wind, so material from the supergiant can be fed directly 
to the compact object through Bondi accretion, or it can pass to 
the compact object via an accretion disk. The latter mechanism 
is typically found in bright X-ray binaries in which the Roche 
lobe overflow of gas from the OB star supercedes the flow of 
accreting matter. For less luminous binaries, the OB star does 
not fill its Roche lobe and the behaviour of the X-ray source 
is determined predominantly by the stellar wind. X-ray emis- 
sion in supergiant HMXB systems is usually persistent, with 
flares stemming from inhomogeneities in the wind. Neutron 
stars with strong magnetic fields develop a hot spot for accre- 
tion which can result in a pulsation. 

The number of persistent and heavily-absorbed HMXBs as- 
sociated with supergiant companions has increased thanks to 
deep, wide-field observations by INTEGRAL combined with 
follow-up X-ray monitoring by XMM-Newton (I Walter et al 
2003 yRodriguez et aljEoollHill et al]l2005l iLutovinov et " 
2005b). So far, these sources hav e been detected p referen- 
tially in the Norma Arm region JWalter et alJ 120041) which 
features high formation rates of OB supergiant stars. It 
is there that INTEGRAL uncovered its first new source, 
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Fig. 1. Intensity map of the Norma Arm in the 20-40 keV band accumulating data from revolutions 30-185 for an effective 
exposure time of 670 ks. The image is in galactic coordinates and represents an area 2° tall by 5° wide. The galactic equator 
(b = 0) bisects the image horizontally. Neighboring sources from the catalog of lBird et alJ d2004h are also shown. 



IGR J 163 18 -4848 llWalter etali fcotB). The Norma Arm also 
harbors IGR J 16393 -4643. This object was initially discov- 
ered in the X-ray band and lis ted as AX J 163904- 4642 in 
the ASCA Faint Source Catalog (Sugizaki et al.ll200ll) . Its flat 
power-law slope, its absorption (A^h = 13*2 x 10 22 cirT 2 ), its 
lack of radio emission, and its position in the galactic plane 
compelled the authors to classify it as a HMXB. A non-thermal 
radio counterpart was recently detected in the ASCA error box 
which suggests a dust-e nshrouded microquasar interpretation 
for the HMXB ( Co mbi et al.ll2004 . and would help explain 
the possible associatio n with the unid entified EGR ET source 
3EGJ 1639-4702 faartman et alll999h noted bv lMalizia etall 
(2004). 

We detected IGR J 16393 -4643 in the ISGRI GPS data 
of the Norma Arm's first visibility period, and we ob- 
tained a follow-up observation with XMM-Newton. The set of 
INTEGRAL undXMM-Newton data are introduced in Section|2] 
In Section^ the refined X-ray position is used to locate the 
most likely counterpart at other wavelengths. Timing and spec- 
tral analyses are presented in Sections|4] and [5] respectively. 
Finally, we discuss the nature of the source and we offer our 
conclusions in Section|6] 

2. Observations and Data Sets 

2. 1. INTEGRAL Data and Imaging 

The INTEGRAL data consist of all CP pointings during revolu- 
tions 30-260, as well as pointings that were public by January 

3, 2005, which had the source within the ISGRI field of view 
(FOV). To improve the quality of the output mosaic, pointings 
with exposure times below 1 ks were ignored. The resulting 
data set groups roughly 1500 pointings with an average expo- 
sure time of 2 ks each. 

Version 4.2 of the INTEGRAL Offline Scientific Analysis 
(OSA) software was used to reduce raw data into images. 
Source extraction employe d version 18 of the INTEGRAL 
General Reference Catalog (Ebi sawa et al.lEo03l) selected for 
objects detected by ISGRI. These tools are available to the 



public through the INT EGRAL Science Data Centre (ISDC: 
ICourvoisieretal.l2003h. 

Intensity, significance, variance, and exposure mosaic im- 
ages were constructed from background-subtracted images of 
individual pointings. Figure^ provides an example of an in- 
tensity map of IGR J 16393 -4643 and its vicinity in the 20- 
40 keV band from the pointings of revolutions 30 to 185. 
The effective exposure time is 670 ks after correcting for vi- 
gnetting. Using this image, we extracted a source location of 
R.A. (J2000)= 16 h 39 m 05 s and Dec.= -46°42.3' (26" uncer- 
taint y) wh i ch agr ees with and improves the ISGRI position of 
iBird et alJ ll2004h . The mean flux (20-60 keV) of the source is 
0.73 ± 0.02 counts per second (cps), or 4.9 mCrab, at a signifi- 
cance of 36cr. 



2.2. XMM-Newton Data and Imaging 

XMM-Newton Hansen et alJ Hqo3 IStriideret alJ HqqJ 
Tur ner et alJEool observed IGR J 16393 -4643 on March 21, 
2004, from 08:21:15 to 11:41:15 (UT). We used the Science 
Analysis System (SAS) v. 6.1.0 software to analyse the data 
and to extract the EPIC spectrum. The data were screened 
for background variability by excluding time intervals in 
which the count rate above 10 keV was greater than a selected 
threshold (0.33 cps for EPIC/MOS1, 0.25 cps for EPIC/MOS2, 
and 2.4 cps for EPIC/PN). After screening, the effective 
exposure times were around 8.4, 8, and 7 ks for EPIC/MOS1 , 
EPIC/MOS2, and EPIC/PN, respectively. 

IGR J 16393 -4643 is clearly detected in images taken by 
EPIC. Figure^ presents an image of IGR J 16393 -4643 from 
the EPIC/MOS1 camera. Furthermore, the source coincides 
with the updated ISGRI position and error circle from this 
paper. The refined X-ray position averaged from EPIC/MOS 
and EPIC/PN is R.A. (J2000)= 16 h 39 m 05.4 s and Dec.= 
-46°42'12" (4" uncertainty). 
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Fig. 2. EPIC/MOS1 image of IGR J 16393 -4643 in the 0.15- 
12 keV energy range with an effective exposure time of 8.4 ks. 
The refined ISGRI error circle from this paper is superimposed. 




Fig. 3. Image from the 2MASS survey of the region around 
IGR J 16393 -4643 in the K s band. The XMM-Newton error cir- 
cle from this paper does not intersect the error boxes of the 
radio source MOST J1639. 0-4642 and the far infrared source 
IRAS 16353-4636. The only infrared object inside the XMM- 
Newton error circle is 2MASS J16390535-4642137. 



3. Counterparts 

We searched for counterparts at other wavelengths and found 
a single infrared s ource b elonging to the Two Microns All-Sky 
Survey (2MASS: ICutri et alJ[2003l) in the XMM-Newton error 
box. This potential counterpart, 2MASS J16390535-4642137, 
is located about 2" away from the XMM-Newton position at 
R.A. (J2000)= 16 h 39 m 05.36 s andDec.= -46°42'13.7" (0.06" 
uncertainty). It appears in the J, H, and K s bands with magni- 
tudes of 14.63+0.06, 13.32±0.04, and 12.78+0.04, respectively 
(95% confidence). 

Figure[3] shows that the XMM-Newton position and error 
radius for IGR J 16393 -4643 do not intersect the error boxes 
of the radio source MOST J1639.0-4642 and the far infrared 
source IRAS 1635 3-4636, whic h are the possible counterparts 
proposed bv lCombi et alJj2004 . The tight XMM-Newton error 
circl e excludes optical counterparts from the USNO-B1 cata- 
log ( Mon etet all20 03). radio sources from the Viz ier database 
JOchsenbein et alJ l20o6T) . and ROSAT sources JVoges et alJ 
Il999fc . 

A possible association of IG R J 16393 -4643 with the 
EGRET source 3EG J 1639-4702 faartman et alJ fl999t) has 
been proposed by Malizia et al. ( 2004). Another ISGRI source, 
IGRJ16358-4726, and 4U 1630-47 are just outside the 
EGRET 95% error contour. In the 2° x 5° degree section of 
the Norma Arm presented in Fig.^ the probability to observe 
a counterpart in the 0.56° error radius of the EGRET source is 
close to 1, so the coincidence with IGR J 16393 -4643 is prob- 
ably a chance one. 

4. Timing Analysis 

4.1. Long-term Variability 

Most ISGRI pointings have exposure times that are too brief 
for a significant detection. Therefore, a mosaic was prepared 
for each 3-day spacecraft revolution (from 30-244) in the 20- 



60 and 60-150 keV bands. We extracted the source flux, error, 
and significance by fitting a Gaussian with a fixed point spread 
function to the mosaic images. Upper limits (3cr), calculated 
from the variance maps, are provided when the source is not 
detected at the Act level in a mosaic image with an effective 
exposure time above 7 ks. There were no detections in mosaic 
images of the 60-150 keV band so its light curve is omitted. 

Figure|4] illustrates the persistence of IGR J 16393 -4643 in 
the 20-60 keV energy range. Table^collects the 15 revolutions 
in which the source is detected in the mosaic image with a sig- 
nificance above Act. The average flux in these mosaics is 0.86 
cps (5.6 mCrab) with a cumulative exposure time of 1.5 Ms. 
The mean flux per revolution varies by a factor of at least 6 
from 0.21+0.05 counts per second (cps) or 1.4 mCrab in revo- 
lution 100 (205 ks), to 1.39±0.14cps (9.3 mCrab) in revolution 
55 (96 ks). When the flux value during revolution 100 is com- 
pared to the highest flux registered in a single 2-ks pointing of 
the 20-60 keV band (4.3 ±0.6 cps, or 29 mCrab, during MJD 
52673.623-52673.641), we find that the source flux varies by a 
factor larger than 20. 

4.2. Pulsations 

During the EPIC/PN observation of 8 ks, the source count 
rate varied from 1 to 4 cps (Fig.[3J. Furthermore, the varia- 
tions occur periodically. By searching for modulations in the 
power spectrum and in the x 2 distribution, we obtain a period 
of 912+5 s at ax 2 of -600 for 30 bins per period (Fig.EJi). The 
pulse period is determined by the centre of a Gaussian fit to 
the x 1 distribution, and the error is ca lculated using the proce- 
dure developed by Home & Baliunas ( 1986) on Lomb-Scargle 
periodograms. The EPIC/PN pulse profile folded with a pe- 
riod of 912 s, illustrated in Fig.^b, indicates a pulse fraction of 
38+5%. 



4 



Bodaghee A. et al.: IGR J16393-4643: a new heavily-obscured X-ray pulsar 




XMH 



52700 52B00 52900 53000 53100 53200 
Time [MJD] 



Fig. 4. ISGRI light curve (20-60 keV) of IGR J 16393 -4643 
where each point represents the flux averaged over 1 spacecraft 
revolution. Upper limits (3cr) are provided when the source is 
not detected for an effective exposure time above 7 ks. The 
XMM-Newton observation (Fig.|5j occurred on MJD 53085. 

For ISGRI, the pulse is best detected in the 100-s light curve 
of revolution 38 in the 15^-0 keV band since this revolution has 
enough effective exposure time (50 ks) with low off-axis angles 
(~ 4°), and its light curve is free of gaps that hinder a period- 
icity search. Figure^ presents the x 1 distribution centered at 
912 s, with 9 bins per period. The phase diagram is shown in 
Fig-EJ), and the pulse fraction is 57+24%. We merged the 100-s 

Table 1. Source flux value (20-60 keV), significance, average 
off-axis angle, and exposure time when IGR J 16393 -4643 is 
detected at the Act level in the mosaic image of a spacecraft rev- 
olution. The median MJD of the visibility period is also given 
in reference to Fig.@] 
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Fig. 5. EPIC/PN light curve for IGR J 16393 -4643 in the 2-10 
keV energy range with a time resolution of 10 s. The observa- 
tion starts at MJD 53085.37516. 

light curves (15-30 keV) from revolutions 37-55, and derived 
a pulse period and error of 912.0+0.1 s. 

With respect to the pulse fractions, the amplitude does 
not change significantly as a function of energy. Both display 
a jagged rise to a peak flux followed by a drop. There ap- 
pear to be primary and secondary peaks in the pulse profiles 
of both instruments at phases of ~0.7 and ~0.2, respectively. 
Although the ISGRI and EPIC periods are derived from obser- 
vations about 400 days apart, the period is not accurate enough 
to search for possible variations. In addition, the ISGRI light 
curves of revolutions that are nearly concurrent with the EPIC 
observation present low source significances which makes it 
difficult to search for a period. 

There is an indication in the EPIC light curve (Fig.[5} that 
the source varies over timescales longer than the pulse period. 
At the beginning of the observation, the maximum of the pulse 
is about 3 cps. This maximum rises to 4 cps about 5 ks after the 
start time, and then decreases. 

We did not detect an orbital period of the order of a few 
days in the combined ISGRI light curves (15-30 keV, 100-s 
resolution) of detected revolutions between 37 and 55. 

5. Spectral Analysis 

5. 1. Average Spectrum 

To extract an ISGRI spectrum, we created mosaics accumulat- 
ing the data from revolutions 37-185 in the following energy 
bands; 22.09-25.92, 25.92-30.23, 30.23-40.28, 40.28-51.29, 
and 51.29-63.26 keV. These boundaries were chosen to con- 
form with the current response matrices. The source is clearly 
detected up to 40 keV, so spectral bins at higher energies are 
ignored. A power law fit to the ISGRI spectrum returns a rela- 
tively steep photon index of 4.5+0.4 (xl = 0.32). 

Spectral extraction for EPIC/PN, EPIC/MOS1, and 
EPIC/MOS2 relied on single and double events within a cir- 



Bodaghee A. et al.: IGR J16393-4643: a new heavily-obscured X-ray pulsar 



5 





800 1000 
Period [s] 




Fig. 6. Top (a): Period search (x 2 distribution) on the ISGRI Fig. 7. Top (a): Distribution of x 2 for the EPIC/PN 2-10 keV 

light curve with 30 bins per period, a resolution of 1 s, and 
centred at 912 s (vertical line). Bottom (b): Folded EPIC/PN 
light curve for a pulse period of 912 s beginning at MJD 
53085.38738. The intervals indicate the pulse and plateau 
states of the phase-resolved spectrum. 



light curve (15-40 keV) of revolution 38, centred at 912 s 
(vertical line), with 9 bins per period, and a resolution of 0.4 
s. Bottom (b): Pulse profile of the folded ISGRI light curve 
(15-40 keV) for a period of 912 s. The start time is MJD 
52674.53833. 



cle of radius 25" centred on the source. To estimate the back- 
ground, we made an event list from a circle of equal radius in 
the same detector and at an equivalent distance from the detec- 
tor's edge. Channels were configured to collect 30 counts per 
bin. An absorbed power law fit to the EPIC spectra has a photon 
index of 1.0+0.1 with a column density of 3.1+0.1xl0 23 cm~ 2 
Or* = 1.15). 

The EPIC spectrum is heavily-absorbed below 4 keV. Iron 
emission lines appear at 6.4 keV (FeK ff ) and at 7.1 keV 
(FeKg). A discontinuity in the continuum above 7.1 keV sug- 
gests an iron absorption edge. There is an indication of an ex- 
cess of soft emission between 0.5 and 2 keV. Including a soft 
blackbody component in an absorbed power law adjusted to 
the EPIC spectra decreases the x 2 by 9 for 680 degrees of 



freedom (dof), or a 0.5% probability that this feature is due to 
chance. Using a partial covering absorption instead of a black- 
body raises the x 2 by 60 for the same dof (679). A partially- 
ionised absorber is also insufficient to model the soft excess. 

Simply fitting an absorbed power law to the spectra from 
EPIC and ISGRI yields N H ~ 4 x 10 23 cirT 2 and T ~ 2.3, but 
the fit is poor (xl = 2.1) with significant residuals in the hard 
X-rays. The EPIC photon index of 1 .0+0. 1 is smaller than the 
one observed for ISGRI (4.5+0.4) which indicates a spectral 
break. Such spectral shapes — a flat power law at low energies 
with a high-energy cutoff between 10-20 keV be yond w hich 
the s lope steepens — are typical of X-ray pulsars i White et alJ 
1 19951) . 
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We fit the combined EPIC and ISGRI spectra (Fig. [8} with a 
broken power law (BP), a power law with an exponential cutoff 
energy (CP), and a Compton emission (CE) model (comptt 
within Xspec). The models are modified by a galactic absorp- 
tion in the direction of th e source (fixed at 2.2 x 10 22 cm~ 2 
JPickev & Lockmar]ll 9901) 1. and an intrinsic photoelectric ab- 
sorption with free iron abundances (Zp e ). Soft excess emis- 
sion is represented by a blackbody (with a fixed temperature 
of kTb = 0.06 keV) that is affected by the galactic absorption 
only. Two narrow Gaussians (with widths fixed at 0) describe 
the iron lines. A constant (C/) accounts for the asynchronous 
observations and for cross-calibration uncertainties. 

The BP model (y 2 /dof = 0.95/675) has photon indices of 
Fj = 0.9+^ and T 2 = 4.6+° ^. However, the BP model does 
not constrain the break energy (£break >17keV) nor C/ be- 
cause it requires one more parameter than the other models. 
The CP model Ct 2 /dof = 1.07/676) gives F = 0.8 ± 0.2, 
and a cutoff energy of Zw = 10 ± IkeV, but residuals re- 
mai n at high energy becau se of the low cutoff temperature. 
In ILutovinov et all ( 2005a). an absorbed cutoff power law fit 
to the combined ASCA and ISGRI spectra had a comparable 
cutoff energy of 11 + 1 keV, but the photon index was poorly 
constrained (r = 1.3 + 1.0). 

Parameters from the CE model (yj/dof = 0.95/676) are 
listed in Table|2] The column density (A^h) is (25 ± 2) X 
10 22 cirT 2 . The comptonising medium has an electron temper- 
ature (kT e ) of 4.4±0.3 keV with an optical depth (r) of 9+1. 
The FeK„ line is at 6.41 ±0.03 and has an equivalent width 
of 60+30 eV when measured with respect to the unabsorbed 
continuum, while the energy of the FeKg line is 7.1+0.1 keV 
with an equivalent width of <120 eV. The detection of the 
line is marginal given the uncertainties and its F-test probabil- 
ity of 3%. However, the ratio of iron intensities (F^/Fko) is 
consiste nt with the value exp ected from the photoionisation of 
iron JKaastra & Mewel|l993l) . The absorbed, integrated fluxes 
(in units of 10 _11 ergscm _2 s _1 ) are 4.4 in the 2-10 keV band, 
and 5.1 in the 20-60 keV band. The observed soft excess flux 
between 0.5 and 2 keV is 7 x 10~ 15 ergscirT 2 s _1 . Since Cj ~ 1, 
the EPIC observation occurred during a period in which the 
source was in an average state. 

5.2. Phase-resolved Spectrum 

The EPIC/PN spectra were binned according to two states: a 
pulse state around the primary peak in the folded light curve 
(phase 0.43-0.77 in Fig. 05), and a plateau state (phases 0.00- 
0.07 and 0.87-1.00). Average count rates for the pulse and 
plateau states are almost a factor of 2 apart at ~1.2 and ~0.7, 
respectively. The CE model was applied to the phase-resolved 
EPIC spectra constrained with the ISGRI spectrum. Bins for 
the plateau gather 20 counts and bins for the pulse collect 40 
counts per channel. The soft excess is not prominent in either 
spectrum so the blackbody component is omitted. This prevents 
an evaluation of the influence of the pulse on the soft excess. 
Table|2l reports parameters for the CE model fit to the phase- 
resolved spectra. The pulsation affects the normalisations but 
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Fig. 8. Photon spectrum of IGR J 16393 -4643 from EPIC/PN, 
EPIC/MOS1, EPIC/MOS2, and ISGRI, fit with the comp- 
tonised continuum model (CE). 
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Fig.9. Phase-resolved spectrum of IGR J16393-4643 from 
EPIC/PN during the pulse and plateau states. The ISGRI spec- 
trum is included to constrain the CE model. 

does not modify the shape of the spectra nor its parameters, 
specifically the Nh, in any appreciable way. 

6. Discussion 

Observations of IGRJ16393-4643 by INTEGRAL and XMM- 
Newton present a source that is highly-obscured and persistent, 
with an intensity varying by a factor larger than 20. The source 
has a pulsation period of 912.0+0.1 s. These attributes, along 
with its spectral shape, its neutral iron lines, and its lack of 
radio emission, suggest a HMXB system consisting of a mag- 
netised neutron star orbiting inside the dense stellar wind of a 
supergiant companion. In the 2-10 keV energy band, the un- 
absorbed flux of IGR J16393-4643 is 9.2 x 10" 11 ergs cirrV 1 . 
Assuming a luminosity of 1.2 x 10 3 6 ergss~ 1 , which is typi- 
cal for accretion-driven X-ray pulsars ( Bild sten et alJ1 997). the 
source is located approximately 10 kpc away. 

The absorbing column density is large whenever the source 
is detected, whether by ASCA (Nh = 13^ x 10 22 cm" 2 ), or now 
with INTEGRAL and XMM-Newton (N H = (25+2)xl0 22 cirT 2 ). 
The strong absorption below 5 keV is intrinsic since it is an or- 
der of magnitude larger than the galactic absorption along the 
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Table 2. Parameters from the Compton emission (CE) model fit to the spectra of EPIC and ISGRI combined. Also listed are 
the values of the CE model (without a blackbody) fit to the pulse and plateau states of the phase-resolved EPIC/PN spectra 
constrained with the ISGRI spectrum. The integrated fluxes are listed as observed (F2-10) or unabsorbed (^ A 10 ). Errors denote 
90% confidence. 



Parameter CE Pulse Plateau Unit 



c, 


0.8+0.2 


0.6+0.2 


1.3+0.7 




T 


9±1 


12+6 


8+ 4 






4.4±0.3 


4.4+0.3 


4.4+0.4 


keV 


N H 


25±2 


23 +2 

^ -6 


24!' 


10 22 cm- 2 


ZFe 


1 +04 


9 +0 - 6 


1 9+0.8 
-0.6 


z G 


^2-10 


4.4 


5.2 


3.0 


10" 11 ergs cm -2 s" 


rtJA 
^2-10 


9.2 


9.8 


6.7 


10~"ergscm~ 2 s" 


Eato 


6.41±0.03 


6.44+0.03 


6.41+0.04 


keV 


FKa 


1.1+0.2 


1.1+0.4 


0.8+0.4 


10- 4 phcnr 2 s" 1 




7.1+0.1 


7.1+0.3 


7.2+0.3 


keV 


Fkp 


<0.62 


<1.1 


<1.3 


lO^phcnT 2 s -1 


Xv/d-o.f. 


0.95/676 


0.92/140 


0.92/94 





line of sight (Nn = 2.2 x 10 22 crrT 2 ). The absorption does not 
vary with the pulsation, which means that it is not related to 
the accretion column, but it may change with the orbital phase. 
This absorption indicates the presence of optically-thick mate- 
rial surrounding the compact object, which is consistent with 
the detection of iron lines. 

Emission lines at 6.4 and 7.1 keV trace the quantity of mat- 
ter in the shell that envelopes the neutron star. These lines are at 
the positions that would be expected from the fluorescence of 
cold neutral matter illuminated by continuum X-rays, and their 
equivalent widths are compatible with a spherical distribution 
of dense gas around the compact object ( Matt 2002). A Ka line 
at 6.41+0.03 ke V (> 1.925 A ) corresponds to an ionisation of at 
most Fe XVIII (House 1969), which does not constrain the dis- 
tance between the ionising source and the inner surface of the 
ionised shell. Given the errors on the line equivalent widths, 
we can not determine whether the fluxes of the iron lines are 
modified by the pulse. 

There appears to be a soft excess flux of 7 x 
10~ 15 ergscm~ 2 s _I that is best represented by a blackbody 
rather than with a partial covering or partially-ionised absorber. 
The scattering of X-rays by the stellar wind i s the most li kely 
expla nation for the soft excess emission jHaberl & W hite 
Il990l) . To what extent the pulse affects the soft excess emis- 
sion is still unknown as the soft excess flux is not prominent in 
the phase-resolved spectrum. 

The ASCA Faint Source Catalog dSudzaki et alJ l200ll) 
listed IGR J 16393 -4643 among the brightest objects in its cat- 
alog, but the uncertain ty of the positio n and the heavy ab- 
sorption prevented Sugizaki et al. (2001) from associating this 
source with an optical counterpart. Malizia et alJ (120041) noted 
the possible association between IGR J 16393 -4643 and the 
unide ntified EGRET source 3EG J1639-4702 faartman et alJ 
Il999i) . Non-thermal radio emission was recen tly detected 
within the ASCA error box dCombi etalJl2004 which sug- 
gests a dust-enshrouded microquasar and could help justify 
the IBIS and EGRET association. However, the XMM-Newton 



position (R.A. (J2000)= 16 h 39 m 05.4 s , Dec.= -46°42'12", 
4" uncertainty) is incompatible with the position of the ra- 
dio source MOST J1639. 0-4642, or any known radio source 
from the Vizier database. The pulsation and the lack of ra- 
dio emission invalidate the microquasar interpretation for the 
HMXB. Also, the chance association of the EGRET source 
with IGR J 16393 -4643 is strong given the high density of 
sources in the region. 

A single infrared candidate, 2MASS J16390535-4642137, 
lies within the XMM-Newton error circle of IGR J 1 6393 -4643 . 
The high-energy spectral and timing characteristics of 
IGR J 16393 -4643 are reminiscent of other heavily-absorbed, 
wind-accreting X-r ay pulsars with OB supergiant companions 
llWalter et al]l2004T) . A Be stellar companion is unlikely given 
that such systems are usually transient. Infrared observations 
should confirm the supergiant nature of the companion. If this 
is the case, it will constitute another argume nt rejecting the 
EGRET association (Orellana & Romero 2005). 

Recently, combined INTEGRAL and XMM-Newton ob- 
servations of the Norma Arm have revealed more pulsat- 
ing X-ray binaries than were previously known. These ob- 
jects have joined the ranks of what might be a new class of 
heavily-absorbed HMXBs that were previously undetected be- 
low 5 keV. The increasing sample size these objects represent 
should enable meaningful statistical studies to be performed. 
Understanding the nature of sources such as IGR J16393-4643 
could shed light on the structure of stellar winds, provide con- 
straints on the masses of neutron stars, and help elucidate the 
evolution of binaries. 
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